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Role of inter-edge tunneling in localizing Majorana zero modes at the ends of quasi 

one-dimensional p + ip systems 

Jimmy A. Hutasoit 1 ^ and Ajit C. Balram 1 '^ 

department of Physics, The Pennsylvania State University, University Park, Pennsylvania 16802 

(Dated: April 11, 2013) 

Potter and Lee have demonstrated the presence of Majorana zero modes at the ends of quasi 
one-dimensional (1-D) p + ip superconductors [T|. We use conformal field theory (CFT) methods to 
show that inter-edge tunneling of the vortex excitations along the length of the channel is crucial for 
such localization. We show that localization of Majorana modes occurs also in quasi 1-D channels of 
the 5/2 fractional quantum Hall (FQH) systems when modeled, following Moore and Read, as the 
p-\-ip paired state of composite fermions. We propose a tunnel-interferometry experiment to detect 
these modes, which should show a tv phase shift of oscillations depending on whether a localized 
Majorana zero mode is present or not, which, in turn, can be controlled by varying the tunneling 
strength. 
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Recently, the zero energy modes of Majorana fermions 
have been receiving a substantial amount of interest due 
to their potential applications in decoherence-free quan- 
tum computation [2] (for recent reviews, see for exam- 
ple [SI [4]). The proposals on how one can realize Majo- 
rana zero modes in condensed matter systems include a 
superconductor-strong topological insulator heterostruc- 
ture [5 and semiconductor-superconductor heterostruc- 
tures [6j [7] . The basic idea of these proposals is to pro- 
duce 1- or 2-D p + ip superconductors, wherein the Ma- 
jorana zero modes will be localized either at the ends of 
the wire [2] or in the vortex cores [SJ [9] . Numerical cal- 
culations show that Majorana zero modes at the ends of 
the p + ip superconducting wire are robust beyond the 
strict 1-D limit as long as the sample width is not much 
larger than the superconducting coherence length [1 . 



One aim of this Letter is to show, using CFT meth- 
ods, that the crucial ingredient of localizing Majorana 
zero modes at the ends of quasi 1-D p + ip supercon- 
ductors is the inter-edge tunneling, which is dominated, 
in a renormalization group (RG) sense, by the tunneling 
of the vortex modes. Since the edge theory of a p + ip 
superconductor also describes, according to the Moore- 
Read model, the neutral sector of the v = 5/2 FQH sys- 
tem, we will use analogous CFT methods to show that 
Majorana zero modes also exist at the ends of quasi 1-D 
v = 5/2 FQH systems. Finally, we propose that a geome- 
try with variable width may be able to detect such modes 
through a tt phase shift in an interference experiment. 



can tunnel across this constriction 1 . Such a constriction 
can be obtained by directly manufacturing a long strip 
with small width or by applying a gate voltage on ei- 
ther sides of the superconductor [10] . Allowing the edge 
excitations to tunnel across the superfluid will be seen 
as crucial in producing Majorana zero modes at the end 
of the constriction. Inter-edge tunneling will not only 
open a gap at the constriction, but also cause the lowest 
energy configuration to contain vortices at the ends, im- 
plying localized Majorana zero modes at the ends of the 
constriction. 

The edge theory of a p+ip superconductor is described 
by a deformed chiral Ising CFT and, focusing on the 
constriction part, the (Euclidean) action is given by 
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= dxdr (£ L + C R + Aim) , 
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where £l and £r refer to left and right moving Majorana 
modes, the "deforming" term C tun describes the tunnel- 
ing, and r = it is the imaginary time. The conformal di- 
mensions of the fields are [iPl,r] = 1/2 and [ctl,r\ = 1/16. 
Here, cjl r are chiral spin fields 2 . 



We start by considering the p-\-ip superconductor. We 
find it convenient to consider a dumbbell-like geometry 
depicted in Fig. [I] The key idea is to have a long constric- 
tion whose width is not much larger than the supercon- 
ducting coherence length, such that the edge excitations 



This is a necessary but not sufficient condition for such tunneling 
to occur, as the presence of disorder is also needed. 
Loosely speaking, the chiral spin field, which corresponds to a 
vortex, can be thought of as the "square root" of the non-chiral 
spin field. This statement is not exact due to the fact that the 
non-chiral spin field does not factorize into the product of chiral 
spin fields |10| . 




FIG. 1: The dumbbell-like geometry of the p + ip su- 
perconductor. The solid line represents the chiral edge 
states and the dashed lines represent the tunneling of the 
edge excitations across the constriction. 

The tunneling operators in Eq. [4] are relevant in the 
RG sense and the RG flow is dictated by 
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dl 



dl 



(5) 



These RG equations are not valid up to arbitrary strong 
coupling as one expects the "weak" tunneling of Eq. [4] 
will become unstable and will result in breaking the su- 
perfluid into multiple weakly connected droplets [10] 

In this regime of "weak" tunneling, the tunneling term 
involving spin fields is more relevant than the one involv- 
ing Majorana fields and thus, A^ <C A a . Furthermore, 
from the first RG equation, we see that the less relevant 
term, which can also be thought of as a mass term, re- 
sults in opening a gap along the constriction because A^ 
grows as one flows to the infrared. 

To understand the effect of the more relevant term, it is 
instructive to continuously deform the (1+1)-D manifold 
of the edge theory into the geometry depicted on the 
first line of Fig. [5J Using the sewing property of (1+1)- 
D CFT, we obtain the second line of Fig. [2j where we 
have factorized the manifold by inserting complete sets of 
states at the cuts 3 . The middle part, which corresponds 
to the constriction, is nothing but the transfer matrix 
element along the strip, which in diagonal basis is given 
by Z^ d ^> = 5y d y exp (— ZJI^ d ), where L is the length 
of the constriction and the spectrum H^ d is obtained by 
diagonalizing the Hamiltonian [T2HT5] . If L ^> 1 is the 
largest length scale in the system, then the sum over the 
states is dominated by the highest weight state \Q) of the 
deformed CFT. 

When there is no tunneling, the edge theory at the 
constriction is just the usual Ising CFT and thus, |0) = 
|0), where |0) is its vacuum state. However, as we will see 
below, tunneling causes \Q) to be a linear combination of 
Ising CFT vacuum state |0) and the so-called spin state 
|<j), which contains a vortex. 

To obtain |fi), we need to calculate the spectrum of the 
deformed CFT and find the state with the lowest energy. 
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FIG. 2: The edge theory lives in a (1+1)-D spacetime 

that looks like a strip sandwiched by two cylinders. 

Here, the vertical direction is the temporal direction. 

Using sewing property of (1+1)-D CFT, we factorize 

the geometry by inserting complete sets of states at the 
cuts, which are represented by the dashed lines. The 
middle part is the transfer matrix element along the 
strip Zw and for long constriction, the sum over all 

states is dominated by the state |0) whose energy is the 
lowest in the deformed CFT. 



The spectrum is obtained by diagonalizing 4 

<tt|iJ|*'> = (tt|(L + Lo-^)|tt')+M*l™l*'>. 

(6) 

where \I/ and ty' are the states of the undeformed CFT, 
Lq and L are the Virasoro generators, and c is the cen- 
tral charge, which for the Ising CFT equals 1/2. Here, 
following Ref. [12], we have kept only the most relevant 
term that deforms the theory away from the conformal 
fixed point. Due to the integr ability of a certain class 
of unitary CFTs, which includes Ising CFT, the above 
statement is not a perturbative statement [12] and it is 
valid even for large X a . Instead, the accuracy of the spec- 
trum obtained depends on how many states one includes 
before diagonalizing the matrix (ty\H\^ f ) [T3HT5] . How- 
ever, since we are interested in knowing only what the 
state with the lowest energy is, it suffices to include only 
the primary states of the Ising CFT: |0), \a) and |e). Sim- 
ilarly, the small corrections due to the less relevant tun- 
neling term will not affect our qualitative results. Here, 
\a) and |e) are defined as per usual radial quantization 
procedure [16] as 

W) = a nc (z = 0,z = 0)|0), |e) = e(z = 0, z = 0)|0), 

(7) 



3 Such property has been used to calculate higher order amplitudes 
in perturbative string theory, see for example Ch. 9 of Ref. |11| . 



4 We assume periodic boundary condition along the temporal di- 
rection T ^ T + 27T. 



where a nc (z,z) and c(z,z) are the non-chiral spin field 
and the energy field, respectively. We have also intro- 
duced complex coordinate z = exp(— x — ivr) and its 
complex conjugate z. The energy field is the product of 
the chiral Majorana fields e(z,z) = ii/j(z)ilj(z), while the 
non-chiral spin field does not factorize into the product 
of the chiral spin fields [10]. Here, we have also replaced 
the (left or right) chiral notation for the fields with holo- 
morphic or anti-holomorphic notation. 

The non-vanishing off-diagonal matrix elements are 
then given by 

(0|<xa|<x) = (a(l)a(l)a nc (0, 0)), 

{a\aa\e) = lim z'lz'l {a nc {l/ z, l/z)a(l)a(l)e(0,0)), 

z,z^>0 

(8) 

and their complex conjugates. We can calculate these 
2- and 3- point functions by introducing another copy of 
the theory, bosonizing and taking the "square root" of 
the corresponding boson correlators [T7] . The non-chiral 
spin field is bosonized such that [17] 

(.-.(j nc (z,z)-.-) 2 - (-■■(j* c (z,z).-.)(.-.(j b nc (z,z)--.) 

(j)(z,z) 



= N 



• cos 



(9) 



where a, b label the different copies of the theory, 
<j>(z, z) = cj)(z) - <j>{z), {<j>(z, z)<j>(w, w)) = - log \z - w\ 2 , 
and the normalization constant N is chosen such that 
the correlation functions will exhibit the correct short 
distance behavior. For example, the 2-point function is 
given by 

<t>{z,z) <j)(w,w) 
— - — cos 



(cr nc (>, z)cr nc (w, w)) = 2 ( COS 



(10) 



Similarly, for the tunneling operator, we have 

(• • • a(z)a(z) • • • ) 2 = N(- • • e ±l ^ z ^> 2 • • • ), (11) 

where the normalization constant N is again chosen such 
that the correlation functions will exhibit the correct 
short distance behavior, while the relative signs between 
tunneling operators in the correlation functions are cho- 
sen to match the chiral conformal blocks [10] . For exam- 
ple 

\2 



{a{z)a{z)(i{w)a{w)Y = ( e ^^ e -^^/^ ( 12 ) 



(z,z)/2 p-icf>(w, w)/2 

in agreement with Refs. [10J[T8] 5 . We then have 



(a(z)a(z)a nc (w,w)) 2 - N ( e^*-*)/ 2 cos ^'^ 



(13) 



We can evaluate TV from the operator product expansion 
(OPE) of a(z 1 )a(z 1 )a(z 2 )a(z 2 )a nc (z 3 , z 3 )a nc (z 4 , z A ) 
which can be expressed as either 



cr(zi)cr(zi)cr nc (z 3 ,Z 3 ) 

a(z 1 )a(zi)a(z 2 )a(z 2 ) 
find N = 2. 

Therefore, we have 



o-(z 2 )o-(z 2 )cr nc (z4, z 4 ) 

0-nc(z3,Z 3 )a nc (z^Z 4 ) 



or 
We 



(0\aa\a) = (a(l)a(l)a nc (0, 0)) 



{(r\aa\e) = lim z~Sz~s(a nc (l/z,l/z)a(i)a(l)e(0,0)) 

z,z^>0 
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In the basis {|0), |<r), |e)}, we have 

(*\H\V'} = ( X, 



(14) 
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and 



|fi) = \0)-8\ a \a) + O(\ 2 a ) 



(15) 



(16) 



Since |0) now contains |cr), the theories living on the 
cylinders (see Fig. [2| are chiral Ising CFTs, but with 
spin fields inserted at the cuts. The effect of the spin 
field on the Majorana field can be seen from the OPE 



^(w)a nc (z,z) 



3 i7r/4 



V2(w-zy/ 2 



— \i{z,z\ (17) 



where the disorder field \i has the same dimension of 
a nc . From the above equation, we see that if we ro- 
tate z around w by an angle of 2n we pick a phase of 
[exp(±27rz)] 1//2 = — 1. Thus, the spin field introduces a 
square root branch cut in the fermion correlators. In 
other words, inserting the spin field <j nc (z,z) at a given 
point in spacetime changes the fermionic boundary con- 
ditions around the point from antiperiodic to periodic, 
and thereby introduces a Majorana zero mode that is lo- 
calized at that point. Therefore, the long constriction 
exhibits localized Majorana zero modes at its ends, in 
agreement with Ref. pQ. 

In the case of v = 5/2 FQH systems, the edge the- 
ory is described by a CFT which consists of chiral Ising 
CFT and a scalar field cp. Again, focusing only on the 
constriction part, the (Euclidean) action is given by 



S = dxdr (jC l + £ R + £ tun ) , 



(18) 



5 Refs. |10l I18| give the chiral correlation functions, from 
which one can deduce the anti-chiral correlation func- 
tions. In particular, {a(z)a(z)a(w)a(w)) 2 is matched with 
(a a (z)a b (z)a a (w)a b (w)) (a a (z)a b (z)a a (w)a b (w)) of Ref. [To]. 
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The RG flow of this theory is given by ([exp(iipL r)} 

1/2) 



dl 
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As is in the p + ip superconductor case, the tunneling 
term involving chiral spin fields is the most relevant term. 
Moreover, from the first two RG equations, we see that 
the other tunneling terms result in opening gaps for the 
edge modes, both the charged scalar modes and the Ma- 
jorana modes. Just like in the p+ip superconductor case, 
the most relevant tunneling term in this case will also re- 
sult in the insertion of non-chiral spin operators at the 
end parts of the geometry. 

The primary states of the scalar sector of the theory 
that are relevant to our problem are given by 



± cp) =: exp I db i 
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where : : denote the normal ordering. In the basis {| 
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and 



\Q) = \0)-8X a \a)+O(Xl) 



(24) 
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Summary and Discussions - Combining the results of 
the RG equations and the deformed CFT spectrum on 
a strip, we have shown that the tunneling of the edge 
excitations across the p + ip superfluid or v = 5/2 FQH 
liquid in the constriction results in localized Major ana 
zero modes at the end parts of the geometry. We can ob- 
serve these zero modes using an interferometry measure- 
ment as depicted in Fig. [3] Here, we assume a "sharp" 
edge, where the electrons from the leads can tunnel all 
the way into the half-filled first excited Landau level. Let 
us consider a slab of v — 5/2 FQH system with significant 
width and no applied gate voltage Vq — 0, so that there 
is no edge excitations tunneling. We then increase Vq to 
create the constriction and thus, at a certain point, al- 
lowing edge excitations to tunnel across the constriction. 
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FIG. 3: As one increases Vg, the FQH liquid is pinched 

to create the constriction (depicted as transitioning 

from the lighter shade to darker ones) and thus turning 

on Majorana zero modes at the ends. 



This results in having Majorana zero modes at the ends, 
which introduce a tt phase shift in the interferometric 
oscillation measured through the leads at the end part. 

As we increase Vq further, we will eventually cross 
into the "strong" tunneling regime, where we pinch off 
the FQH liquid totally, splitting it into two decoupled 
FQH systems. At this point, there will be no chiral spin 
fields tunneling across [10] and therefore, there will be 
no Majorana zero modes at the ends. This results in a 
phase shift back to the original interferometric pattern. 
To summarize, starting at Vg = 0, as we increase Vg, 
there will be a phase shift in the interferometric pattern 
at some point, which will then shift back as we keep on 
increasing Vg- 

Lastly, let us note that we have treated the v = 5/2 
FQH state assuming the Moore- Read Pfaffian model, al- 
though the presence of Majorana edge modes from inter- 
acting microscopic theory is not yet conclusive [19 . An 
observation of localized Majorana modes will serve as a 
confirmation of certain non-trivial aspects of this model. 
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